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1. Introduction 

Stars with masses in the 1-8 Mq range are very hkely to go through the evolutionary phase of 
PN ejection at the end of their hves. The ejection of most of the gaseous envelope in Asymptotic 
Giant Branch (AGB) stars occurs at the tip of the Thermal Pulse phase, in the form of a massive, 
low velocity wind (Iben &; Renzini 1983). Subsequently, the remnant CS ionizes the gaseous ejecta, 
while the fast, low mass-loss rate CS wind shapes the PN. The observed PN morphology thus 
depends on a complicate combination of phenomena, some occurring within the nebular gas, which 
evolves in dynamic time scale, others deriving from the evolution of the stellar progenitors and of 
the CSs. Morphology may also depend on the physical status of the interstellar environment of the 
PN progenitor. 

The relations between PN morphology and CS evolution may be used to explore the gas 
versus star interplay. Pioneering this field, Peimbert (1978) found a link between PN morphology 
and stellar evolution through the chemical elements that are produced by PN progenitors, and 
are drawn to the stellar surface via dredge-up episodes. Peimbcrt's study showed that most of 
the helium- and nitrogen-rich, carbon-poor PNs where of asymmetric (bipolar etc.) shapes. This 
result is broadly consistent with the predictions of stellar evolution if the progenitors of asymmetric 
PNs have on average larger masses than the progenitors of round (and elliptical) PNs, without any 
assumption or relation to stellar multiplicity of the stellar progenitors. 

About a decade ago, several authors have tried to tie up the correlations between nebular 
morphology and stellar properties (Stanghellini, Corradi, & Schwarz 1993; Corradi Sz Schwarz 
1995), inspired by the new bounty of narrowband images available for Galactic PNs (Schwarz, 
Corradi, &: Mclnick 1992; Balick 1987). Stanghellini, Corradi, &: Schwarz (1993), based on a 
large, yet non homogeneously selected, sample of Southern PNs from the ESO catalog (Schwarz, 
Corradi, & Melnick 1992), found that the CSs of symmetric and asymmetric PNs have different 
mass distribution, and that the CSs of bipolar PNs are spread over a large mass range. Other 
papers have been written in the recent years (Amnuel 1995; Gorny, Stasinska, &; Tylenda 1997) on 
this subject, none of which have challenged the earlier results. 

Strength to the association of bipolar PNs with high mass progenitors comes from the Galac- 
tic distribution of morphological types. Among others, Schwarz, Corradi, &; Stanghellini (1993), 
Stanghellini (1995), and Manchado et al. (2000) found that bipolar PNs lie closer to the Galactic 
plane than round and elliptical PNs, making the high-mass origin plausible. 

The motivation of the present work is to test the correlations between PN morphology and CS 
evolution with a larger and homogeneous PN database. In 1996 the lAC Catalog of PNs (Manchado 
et al. 1996), and with it a complete and homogeneous PN image data set, became available. 

The aim of this paper is to discuss the relations between morphological type and the physical 
parameters that arc linked to stellar populations and evolution. We will examine principally the 
variations of temperature, luminosity, and mass of the CSs. We will also show the variation across 



-3- 



morphological types of elemental abundances and Galactic distribution. The statistics of the neb- 
ular characteristics variations across morphological types will be published elsewhere (Manchado 
et al. 2002). 

In §2 we discuss the sample selection and the morphological classification. In § 3 we describe 
the determination of the stellar temperature and luminosity. In § 4 we present the results of our 
analysis. A discussion of the significance of our findings for the understanding of the evolution of 
PNs and their CSs is presented in §5. 

2. Morphological classification 

The database that we have used in this paper includes the lAC morphological catalog of 
Northern Galactic PNs (Manchado et al. 1996), supplemented by a selection of PNs from Balick 
(1987) and Schwarz, Corradi, k. Melnick (1992). Therefore, our sample includes all PNs in Acker et 
al. (1992) that are larger than 4 arcsec in maximum size, and that are within the observing limits 
of the lAC and La Palma Observatory telescopes (Manchado et al. 2000). 

We classify the PN morphology primary on the basis of the Ha images, and secondary on the 
[O III] 5007 and the [N ll] 6584 images, when available. A thorough discussion on the extended 
morphological classification will be available in Manchado et al. (2002). Here, we briefly describe 
the morphological classes used to analyze the results on the CSs correlations. 

We use the following, simplified morphological scheme: 

• Round (R) PNs. 

• Elliptical (E) PNs. We distinguish between R and E PNs if an axial difference of at least 5 
% is detected. 

• Bipolar (B) PNs. As commonly defined, bipolarity implies at least one pair of lobes and a 
pinched waist. Quadrupolar PNs belong to this category. 



Together with the three main morphological classes above, we also examine the group prop- 
erties of Bipolar Core (BC) PNs. This morphological class identifies those PNs whose external 
contour does not show obvious bipolar lobes, but whose inner structure is clearly suggestive of a 
projected ring. A ring-like structure was detected above the 20 % flux level for our BC PNs to be 

classified as such. This class was firstly used in Stanghellini et al. (1999) to analyze bipolarity in 
Magellanic Cloud PN images acquired before the first HST servicing mission that restored the tele- 
scope optics. Note that Corny, Stasinska, &: Tylenda (1997) use a similar definition for embedded 
bipolars (BE). 
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A few comments on this simplified classification scheme are in order. First note that Stanghellini, 
Corradi, & Schwarz (1993) do not distinguish between R and E PNs. We now know that it is possi- 
ble to separate these classes, and that they have diff'erent properties. Note that our classification of 
round PN is an objective one, based on the measured axial ratio on the 2D images, without making 
assuptions on the 3D structures (note that Soker (2002) uses a different classsification approach, 
and subjectively defines elliptical those PNs that have structures such as knots or filaments, even 
if they have round contours). 

Second, in this paper we subsume the quadrupolar PNs into the bipolar class. Quadrupolar 
(Q) PNs show two pairs of lobes (Manchado, Stanghellini, k. Guerrero 1996) and their formation 
mechanism has to be similar to that of bipolar PNs. Our database contains but a few Q PNs, and 
studying them as a class apart will dilute their statistical significance anyway. 

Third, we do not classify pointsymmetric (P) PNs as a separate morphological type. As pointed 
out by Manchado ct al. (2000), pointsymmetric structures are present in a variety of main body 
morphologies. For this reason, we classified the pointsymmetric PNs as round, bipolar, or elliptical 
depending on the structure of the main body of the PN. Pointsymmetry is detected in a very small 
fraction of the PNs in our sample, and if they were selected in a separate class the results from 
our analysis would not change statistically. Finally, it is worth noting that the bipolar core PNs 
would have been classified as R or E in the usual or the new morphological schemes (Manchado et 
al. 2002). We want to isolate their properties in this paper, and relate them to those of B PNs. By 
doing so, we moderately decrease the R and E PN data sets. 

3. Temperature and luminosity derivation 

The stellar properties of the selected PNs were determined via the Zanstra analysis (Kaler 
1983). The input parameters to the Zanstra analysis were selected from the literature, with the 
endeavor of using the best possible set of input data, chosen as homogeneously as possible. In Table 
1 we list our reference choices for the input parameters. 

The stellar magnitudes used for the analysis were taken, for the most part, from the paper by 
Tylenda et al (1991). 

The H/3 and He ll (A4686) line intensities were drawn, when available, from Cahn, Kaler, & 
Stanghellini (1992), otherwise from Tylenda et al (1991). Additional fluxes and extinction constant 

references arc listed in Table 1. Extinction constants are from Cahn, Kaler, & Stanghellini (1992) 
or, alternatively, from Stasinska ct al. (1992). If neither of these two references had relevant data, 
we used the additional references listed in table 1. 

Nebular radii were measured by us on the original frames. For elliptical PNs, we compute 
equivalent radii as i?E = [(a x &)/4]^^^, where a and h are the elliptical axes. For bipolar PNs, we 
make the approximation that the lobes are two ellipses, and we define i?E = [(o x ^)/4]^/^, where 
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a is the extension of the waist, and b the extension of the bipolar lobes. 

All distances used in this paper are statistical distances. We have calculated all the statis- 
tical distances based on the newly measured radii, and the Cahn, Kaler, Sz Stanghellini (1992) 
distance scale calibration. The 6 cm fluxes used in this calculation were taken from Cahn, Kaler, &i 
Stanghellini (1992) where available, otherwise from the other references in Table 1. To assure that 
the calibration of the distance scale is independent on the morphological type of the calibrators 
used, we have morphologically classified the calibrators in Cahn, Kaler, &; Stanghellini (1992), and 
verified that all morphological types have similar weight in the calibration. We also verified that 
the calibration based on just bipolar and bipolar core PNs is almost identical to the calibration 
based on round PNs. 

All input data are used in our Zanstra code with their own original error bars, as given from the 
original papers listed above and in Table 1. The Zanstra analysis then produces the propagated 
error bars. The statistical distances do not carry error bars, as it is impossible to quantify the 
magnitude of the errors in the assumptions made to scale the ionized masses with the nebular 
distances (see the discussion by Cahn, Kaler, & Stanghelhni (1992)). 

4. Results 

In the following sections we illustrate the results of our investigation. Among the 255 Northern 
Galactic PNs with morphological classification, Zanstra analysis was performed in all the cases 
where the input parameters were available in the literature, and precisely for 54 round, 96 elliptical, 
36 bipolar, and 19 bipolar core PNs. 

In Table 2 we give the results of the Zanstra analysis. Column (1) gives the PN name; column 
(2) gives the morphological type, following the simplified scheme described in §2. Unless otherwise 
noted within the Table the images for morphological analysis were from Manchado et al. (1996). 
Columns (3) through (6) give the Zanstra temperature (and error) and luminosity (and error), in 
logarithmic format. We always give the He II Zanstra temperature if available. Otherwise, we give 
the H I Zanstra temperature. In the cases in which the stellar magnitude was not available, we 
give the cross-over temperature instead (Kaler & Jacoby 1989). The reader should not use cross- 
over temperatures in Table 2 at face value, since they are calculated assuming that the He II and 
H I temperatures are equal. Cross-over temperatures are upper limits for most nebulae except for 
optically thick PNs, and the magnitudes estimated via cross-over analysis are also upper limits. 

Distances, dimensions, and detailed morphology will be discussed in a future paper. 

In Table 3 and 4 we give the averaged physical diagnostics across morphological types that we 
use in the following sections. Table 3 contains the Galactic distribution and Zanstra temperature 
diagnostics (with the a values in square braketts), while Table 4 gives the elemental abundance 
diagnostics. 
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4.1. Spatial distribution 

The apparent spatial distribution of Galactic PNs with different morphology is shown in Figure 
1, where the different symbols indicate morphological types (open circles are round, asterisks are 
elliptical, filled squares are bipolar, and open squares are bipolar core PNs). It appears from 
this plot that the bipolar and bipolar core PNs are found at lower Galactic latitudes, and their 
distribution do not extend, for the most part, away from the disk. Differently, elliptical and round 
PNs are spread all over the 1-b plane. From column (2) in Table 3 we can see that the round PNs 
have a latitude range much larger than that of bipolar PNs. 

We have run the Kolmogorov-Smirnov (K-S) statistics (Presse et al. 1986) on the Galactic 
latitude distribution to investigate the quantitative differences of the different morphology PNs. 
The K-S statistical test give the maximum absolute value of the difference between two given 
distributions, Dmax- The significance level for the null hypothesis that the data sets are drawn 
from the same distribution (P) is also given by the K-S statistics: lower values of P indicate two 
different distributions. Typically, P < 5 x 10~^ indicates that there is a real possibility that the 
two distributions are different, while P=l indicates identical distributions. Larger data sets give 
better results, but data sets as large as N=20 already give a good indication on the statistics. 

By considering the Galactic latitude distributions of R and B PNs, we found that Djnax=l-8 
and P=0.003. The latitude distributions of round and elliptical PNs are not so sharply dissimilar 
(P=0.19). To enhance the statistical significance of the results we grouped B and BC PNs in one 
class (B-FBC). The K-S probability for R and B-I-BC PNs to have similar b distributions is 0.086, 
while this probabilities is 0.53 when we use the E and B-I-BC samples. Our conclusion is that 
there is a sharp difference between the nature of R and B-|-BC PNs. Much more detail on these 
distribution will be presented in a future paper. 

In Figure 2 we show the PNs distribution of the absolute distance to the Galactic plane, \z\, 
for each morphological type. This parameter is calculated with the statistical distance scale, and 
the sample sizes are different than those of Figure 1, since some PNs did not have the necessary 
parameters in the literature for the derivation of their statistical distance. We see in Figure 2 that 
the distribution of R PNs is the broadest, followed by the E, the BC, and the B PNs. 

We have calculated the average scale distance from the Galactic plane for the different mor- 
phological types (column [3] in Table 3). Round PNs have the largest average E PNs follow, 
then BC PNs , and finally B PNs, closest to the Galactic plane. The average scale height of all PNs 
in the sample is 0.44 kpc. Our calculations basically confirm the early results by Schwarz, Corradi, 
& Stanghcllini (1993). It is worth noting that the average distance per morphological class are 
comparable to one another. 

Our results on the Galactic distribution of morphological types show that morphology could be 
used as indicator of stellar population. For example, the scale height for a young disk population is 
similar to that of bipolar and bipolar core PNs (Scalo 1998). The caveat of this conclusion is that 
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there may be a selection effect that limits the observations of low surface brightness PNs toward 
the Galactic plane. 

In order to asses the importance of the selection effects, we checked the H/3 surface brightness 
distributions of the PNs, calculated as F^/vri?^, where the H/3 flux is corrected for extinction and 
the radii are the equivalent radii calculated as described in §3. We found that surface brightness of 
R and E PNs are in the —17 < log SBp < —11.5 range, while B and BC PN's surface brightness 
range within —15.5 < log SBp < —12.5^. The lack of high surface brightness B and BC PNs 
may simply indicate that they evolve faster than E and R PNs, similarly to what observed in the 
Magellanic Clouds (Stanghellini et al. 1999; Shaw et al. 2001). This is another hint that B and BC 
PNs may originate from higher mass progenitors. 

On the other hand, the lack of very low surface brightness B and BC PNs may indicate 
that there may be a selection effect against low surface brightness B and BC PNs in the Galaxy, 
possibly due to the interstellar extinction which is higher for these PNs, given their location within 
the Galactic plane. The average H/3 extinction for each morphological type is given in column (4) 
or Table 3. Since on average extinction does anti-correlates with the distance from the Galactic 
plane, it roughly indicates that interstellar extinction is predominant, and the internal extinction 
is only marginally important in the overall values. 

4.2. Central stars on the Log L — log T plane, and the inferred mass distribution 

Zanstra temperatures of the CSs were calculated as described by Stanghellini, Corradi, & 
Schwarz (1993). The He ii temperatures are the most reliable, as most PNs are optically thick 
to the He II ionizing photons. H I temperatures are generally reliable for the optically thick PNs. 
When using H i Zanstra temperatures wc should keep in mind that it is difficult to define which 
PN is optically thick and which is thin, especially in the case of bipolar PNs. Thus the only reliable 
results are the few temperatures calculated from the He II line intensities. 

In order to disclose all possible relations between PNs and their CSs, we have taken a con- 
servative approach. We first consider the most reliable Zanstra temperature set, using only He II 
Zanstra temperatures with minimal errors (AlogTefj < 0.1). The averages of column (5) in Table 
3 refer to this set. We do not find marked differences among the average temperatures of stars 
hosted by different morphological types. 

Corradi & Schwarz (1995) found an average effective temperature of 142,000 K for stars hosted 
by bipolar PNs, and 99,000 for those within the elliptical (including round) PNs. We do not disclose 
such a difference in the average temperatures. In Figure 3 we show the location on the log T- log 
L plane of the PN CSs, where AlogTegf < 0.1 and AlogL/ Lq < 0.3 (note that the error bars do not 



This results is obtained with the PNs of Table 2, but we eliminated NGC 7017 from this calculation 
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account for uncertainties in the statistical distances). 

The HR diagram distribution of Figure 3 shows that stars hosted by round nebulae are spread 
on the logT - logL plane differently than stars hosted by bipolar and bipolar core PNs. Central 
stars of elliptical PNs cover more or less the whole plane, up to the 0.8 Mq evolutionary track. We 
derive the distribution of the CS's masses for the different morphologies from Figure 3, and we show 
it in Figure 4. Round nebulae appear to have a peak CS mass around 0.55, then the distribution 
declines for larger masses. Central stars of elliptical PNs have a flatter mass distribution. Bipolar 
and bipolar core PNs seem to have a flat mass distribution (better compare the round and elliptical 
distribution with the sum of bipolar and bipolar core, since these two later classes are so sparsely 
populated) . 

Figure 5 is similar to Figure 3, with the difference that we plot the HR diagram location of CS 
PNs whose temperatures have errors smaller than 0.4 dex, and where we include the H i Zanstra 
temperatures. Figure 6 gives the histograms of the mass distributions relative to the sample of 
Figure 5. Here, the low-mass section of the histogram for elliptical PNs is more prominent than in 
Figure 4, but the overall distributions have the same qualitative form. These results are indicative 
of different CS masses for different PN shapes. The statistical K-S analysis applied to this set of 
temperatures (and the relative luminosities) does not disclose enormous differences between the 
effective temperature and luminosity distributions. 

We have analyzed the ratios of He ll to H l Zanstra temperatures across morphological types, to 
have an indication of their optical thickness. Our aim was to probe the preliminary results of Figures 
1 and 3 in Stanghellini, Corradi, &: Schwarz (1993), by using the better and more homogeneous 
sample of Northern Galactic PNs. Our results (column [6] of Table 3) shows that, although there 
are differences in the average Zanstra ratio for the various morphologies, such differences are milder 
that previously found (Stanghellini, Corradi, & Schwarz 1993). We conflrm that, on average, bipolar 
PNs have lower Zanstra ratio than symmetric PNs, but we can not draw the same conclusion for 
bipolar core PNs based on our small sample (see also the higher sigma value). 

4.3. Hints from elemental abundances 

During the evolution of intermediate- mass stars, the chemical enrichment of the outer region 
of the stars occurs in a series of dredge- up and nuclear burning events (Ibcn &: Renzini 1983), that 
enrich the future PN in different ways, depending on the stellar mass. By measuring the CNO 
elemental abundances we can probe stellar evolution. Carbon depletion and nitrogen enrichment 
depend on the progenitor mass, forming a direct connection between observed progenitor mass (i.e. 
Population) and chemical content. 

The argon and neon abundances in PNs are probes of the original chemical mix at the time of 
the formation of the progenitors. In fact, these elements are unaffected during the evolution of the 
stars in this mass range. 
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A large study of the PN abundances and their correlations to shell morphology is in preparation 
(Manchado ct al. 2002), based on abundances of Galactic PNs newly calculated with a selection 
of the available spectral data in the literature. In Table 4 we show the averages of the derived 
abundances (all abundances are with respect to hydrogen, by number), calculated for each mor- 
phological type. Note that we use only reliable elemental abundances (see Manchado et al. (2002)). 
With each entry, in parenthesis, we report the statistical sample used for that particular diagnostic. 

Among the significant results, we find that helium abundances are, on average, higher in bipolar 
PNs than in round and elliptical PNs. The average helium abundance trace nicely the increasing 
asymmetry, from round to bipolar PNs. If we consider R+E PNs and B+BC PNs, the average 
helium abundance of the two groups differ, quite significantly, at the 25 percent level. 

The oxygen abundance averages do not change as much from one group to the next, with the 
exception that the average oxygen abundance for bipolar PNs is slightly higher than for the other 
groups. 

Bipolar PNs have much higher average nitrogen abundance than the other morphological 
types, confirming the results by Peimbert (1978) and others. Bipolar core PNs are mildly nitrogen 
enriched, based on a rather small PN sample. If we compare the R+E and the B+BC groups, we 
obtain that the average nitrogen abundances in the two groups difi^er at the 60 percent level. This 
is a clear signature that bipolar PNs have massive progenitors. 

Other diagnostics do not show remarkable differences, or are statistically insignificant. We do 
not see the trend found in LMC, where B PNs have higher neon and argon abundances (Stanghellini, 
et al. 2000). 

While the abundances used here are probably the best collection available to date, it is very 
difficult to draw sound conclusions on the basis of sparse spectroscopic data with slits collecting 
variable areas of PNs. This is particularly debilitating for multiple shell or very bipolar PNs, as 
their abundance gradients can be substantial. It is worth remarking that, in addition to the above 
difficulties, the abundances are also contaminated by the ISM material, ingested by the PN during 
their evolution (Villaver, Garcia-Segura, &: Manchado 2002). 

4.4. Evolutionary paths 

As seen in the two preceding sections, there is an indication of mass segregation among the GSs 
hosted by symmetric and asymmetric PNs. This has bearings on the possible evolutionary paths 
that lead to the different morphologies. The formation of asymmetric PNs has been ascribed to 
binary evolution, either through common envelope (CE) or via wide binary or planetary companions 
(Sokcr 1997), as well as to the presence of stellar rotation and magnetic fields (Garcia-Segura, 
Langcr, &; Rozyczka 1999). Soker (1997) selected a large sample of PNs, and, based solely on their 
morphology, determines the likeliness that these morphologies originate (1) from the evolution of 
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a single star, (2) from binary evolution, where the member of the binary system go through the 
Common Envelope (CE) phase, (3) from binary evolution of a pair that avoids the CE phase, or (4) 
from the evolution of a star and a planet companion. We have checked Sokcr's determinations of 
the stellar origins of PNs against our morphological sample. We found, based on our classification 
scheme and on his evolutionary scheme, that (1) no B or BC PNs would originate from single 
stars; (2) no round PNs progenitors would have gone through the CE phase; and (3) no R nor 
E PNs would originate through binary stars, without going through the CE phase. If (1) was 
true, then why the CS masses of B and BC PNs are distributed differently than those of R and 
E PNs? Why B and BC PNs would be found at lower scale heights on the Galactic plane, and 
would they have higher helium and nitrogen content? Clearly, a theoretical scenario that explains 
all the observations is not available at this time. In particular, while the binary star scenario can 
be streched to explain the Galactic distribution of bipolar PNs (Soker 1998), much more stellar 
evolution modeling is needed to compare the chemical yields in the cases of single and binary star 
evolution. 



5. Summary and discussion 

The analysis of Northern Galactic PNs have produced interesting results. First, there is a 
general indication that bipolar and bipolar core PNs belong to the same physical class, even if 
they may be at different evolutionary stages. This results from the apparent and physical spatial 
distribution of B and BC PNs, from the excinction analysis, and from the CSs distributions. The 
BC sample is barely large enough for accurate statistical analysis, thus this result is still to be 
confirmed in the future. It agrees, though, with the analysis of Magellanic Cloud PNs, as discussed 
by Stanghellini et al. (1999); Stanghellini, et al. (2000); Shaw et al. (2001). 

Second, we have analyzed the CSs via Zanstra analysis. The results are in moderate agreement 
with those of Stanghellini, Corradi, & Schwarz (1993), after correcting for the different morpho- 
logical classification. Although the mass histograms of the CSs of Figures 4 and 6 tell us that 
bipolar (and bipolar core) PNs have a larger fraction of massive CSs (i.e., massive progenitors) 
than round and elliptical PNs, the result is not as sharp as previous results have shown. Our 
approach is accurate, yet it lacks true luminosity error bars, as the distances to the nebulae were 
derived in a statistical way. Furthermore, the Zanstra analysis implies the assumption that the CS 
radiates as a blackbody. It is worth noting that, to date, no CS mass has been derived without 
strong assumptions on the distances, on the ionized masses, or on model-dependent parameters. 
Our conclusion is that, to date, there is not a reliable way to measure the mass of Galactic PNs or 
their stellar progenitors, although our analysis seems to indicate that asymmetric PNs were formed 
from more massive progenitors, on average. 

Third, we have investigated the progenitors' populations of the different morphological PN 
types: higher mass progenitors should have different chemistry both from the evolutionary dredge- 
ups and from the different chemistry at star formation. We found that the helium and nitrogen (and 
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very marginally oxygen) average abundances in bipolar and bipolar core PNs are larger than the 
corresponding values in round and elliptical PNs. In agreement with the theory of stellar evolution, 
we confirm that B and BC PNs have more massive progenitors than E and R PNs. 

On the other hand, the average neon and argon abundances are almost identical in roTind, 
elliptical, and bipolar PNs (bipolar core PNs do show high neon and low argon abundances, but 
the sample size is very small for statistical significance). This result is in contrast to a study of 
LMC PNs (Stanghellini, et al. 2000), where bipolar (and bipolar core) PNs have higher neon and 
argon abundances than round and elliptical PNs. We know that there is a larger fraction of bipolar 
PNs in the LMC than in the Galaxy (Shaw et al. 2001), and that Galactic bipolar PNs tend to 
be located in the Galactic plane, where the interstellar extinction is maximum. It seems then 
reasonable that the observed Galactic bipolar PNs are local, so their neon and argon abundances 
reflect the local environment at formation, rather than the average Galactic disk environment. 

Last, we conclude that the study of the statistical correlations between CSs and PN morphology 
has probably reached, in this form, its exploitation. By increasing the statistical significance and 
the quality of the analyzed sample we do not get sharper results on the statistical correlations 
between CSs and PNs. The future of this type of studies is, in our view, two-fold. On one hand, 
the acquisition and interpretation of CS spectra, and at the same time, the development and 
improvement of non-LTE stellar models, is essential to acquire the necessary stellar data for a non- 
biased analysis. On the other hand, better nebular analysis on the individual objects will improve 
the knowledge of reddening to the individual PNs, and their global elemental abundances, possibly 
extending the analysis to UV spectroscopy, i.e., to reliable carbon abundances. High resolution 
spectroscopy on the PNs will also be essential to confirm the morphology of some uncertain cases, 
where the spatial dimension of the velocity is essential. 

A continuing study of those extra- Galactic PNs whose morphologies are detectable with today's 
technology (e.g., the LMC, the SMC, and the nearby dwarf galaxies observed with HST) will have 
the double benefit of minimizing the ill-effects of the unknown Galactic PN distance scale, and 
those of the poorly known reddening, and of allowing us to extend and integrate our knowledge 
of the final phases of intermediate-mass stars to Galaxies of different types and metallicities. The 
emerging picture is that PNs generated by more massive progenitors, those that have been ejected 
by disk population stars, seem to be bipolar in shape, while smaller mass stars give birth to round 
PNs. Elliptical PNs are produced by stars in a wider mass range. This result is in agreement 
with the findings of the PN distribution within the Galaxy, and with the relative N/H and He/H 
abundance of the Galactic PNs sample. It bears on the formation of PNs of different morphological 
types, and it agrees with bipolar PNs to be formed by the evolution of larger mass stars than round 
and elliptical PNs. While both single and binary star evolutionary models can account for bipolar 
PN formation, via common envelope evolution, magnetic fields, and fast stellar rotation, the binary 
scenario alone would not comply with our results. 

It is a pleasure to thank Dick Shaw, Bruce Balick, Stacy Palen, Chris Blades, and Max Mutchler 
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for participating in our collective PN morphology classification experiment. Mark Dickinson is 
thanked for his help in the statistical analysis. 
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Table 2. Morphology, temperatures, and luminosities 



name morph. log Teff Aiogx log L/Lq AiogL 



A4'» 


R 


5.114 


0.302 


1.717 


0.308 


A8* 


R 


5.187 


0.301 


1.932 


0.302 




R 


5.214 


0.301 


2.497 


0.302 


A16 


R 


4.936 


0.040 


1.800 


0.142 


A20'' 


R 


5.012 


0.033 


2.744 


0.105 


A28 


R 


4.967 


0.301 


1.269 


0.305 


A30 


R 


4.874 


0.013 


2.600 


0.068 


A33 


R 


4.979 


0.022 


2.257 


0.063 


A39 


R 


4.953 


0.021 


1.934 


0.078 


ASO^- 


R 


5.178 


0.301 


2.027 


0.301 


A53='''' 


R 


5.215 


0.306 


2.356 


0.333 


A71 


R 


5.113 


0.015 


1.830 


0.042 


A83=' 


R 


4.386 


0.211 


2.082 


3.513 


BV5-3 


R 


4.925 


0.026 


2.447 


0.113 


Bal 


R 


5.058 


0.008 


2.126 


0.025 


CnS-l'' 


R 


4.809 


0.070 


3.893 


0.281 


H3-29^ 


R 


5.954 


0.301 


2.496 


0.301 


H3-75^ 


R 


5.099 


0.301 


1.987 


0.301 


Hel-4'= 


R 


5.133 


0.053 


3.300 


0.106 


Hel-5'' 


R 


4.894 


0.301 


0.000 




HC2-432'' 


R 


4.305 


0.190 


4.713 


3.159 


IC 1454<:'d 


R 


5.067 


0.301 


1.760 


0.304 


IC 4593^^ 


R 


4.450 


0.014 


3.472 


0.041 




R 


5.038 


0.301 


2.100 


0.301 


Kl-9'» 


R 


5.194 


0.301 


1.295 


0.302 


Kl-14 


R 


4.823 


0.021 


3.351 


0.117 


Kl-20* 


R 


5.011 


0.301 


1.349 


0.301 


K3-2=>- 


R 


4.267 


0.178 


4.790 


2.953 


K3-27'' 


R 


5.494 


0.301 


2.499 


0.303 


K3-51'' 


R 


5.534 


0.301 


2.946 


0.303 


K3-56^ 


R 


5.449 


0.301 


0.000 




K3-7'' 


R 


4.179 


0.145 


1.931 


2.365 


K3-73 


R 


5.049 


0.008 


1.779 


0.022 


K3-81 


R 


4.814 


0.012 


3.641 


0.069 


KjPnl^ 


R 


5.156 


0.301 


2.055 


0.301 


LSAia 


R 


4.316 


0.193 


2.926 


3.212 


Ml-77'' 


R 


4.320 


0.006 


4.430 


0.030 


Ml-80* 


R 


5.226 


0.301 


2.748 


0.302 


M4-18'' 


R 


4.432 


0.010 


3.576 


0.034 


NGC 2242 


R 


5.025 


0.033 


2.778 


0.100 


NGC 2438 


R 


5.175 


0.082 


2.055 


0.134 


NGC 3587 


R 


5.049 


0.018 


1.606 


0.053 


NGC 6852 


R 


5.130 


0.014 


2.446 


0.035 


NGC 6879 


R 


4.835 


0.026 


3.617 


0.130 


NGC 6884 


R 


4.902 


0.057 


3.273 


0.226 


NGC 6894a 


R 


4.991 


0.302 


2.444 


0.306 


NGC 7094 


R 


4.874 


0.009 


3.058 


0.042 


NGG6842 


R 


4.994 


0.031 


2.782 


0.105 
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Table 2 — Continued 



name morph. log T^ff ^Io^t log L/L© ^lo^L 



Nal 


R 


4.899 


0.025 


3.597 


0.118 




r> 
K 


4.315 


0.009 


3.113 


0.044 


Sdl^ 


R 


4.301 


0.189 


2.084 


3.137 


k)to-i 


K 


5.134 


0.301 


2.144 


0.301 


vyz-d 


■D 


A '7£iA 

4.7d4 


0.013 


4.153 


0.047 


Wei- 5 


■D 


4.945 


0.011 


2.904 


0.042 


A O 

Az 


17 

Cj 


5.124 


0.030 


2.431 


0.081 




Cj 


5.377 


0.304 


2.306 


0.322 


A 1 Q 
Aio 


17 


c; 1 TA 
O.i / 4 


n 1 1 A 




n 1 01 
U.lol 


A18 


17 

Cj 


5.050 


0.019 


1.818 


0.054 


A19^ 


17 


4.310 


0.191 


1.330 


3.181 


A21^ 


E 


5.062 


0.079 


0.932 


0.045 


A /IQ 

A4o 


17 


4.833 


0.006 


2.778 


0.032 


A4o 


17 

Ej 


4.859 


0.079 


2.169 


0.268 


A49^ 


T7 

Cj 


4.303 


0.190 


2.393 


3.148 


A04 


17 


4.0UO 


n 1 on 


z.Zoo 


o.loy 


A57*^ 


17 
Cj 


4.667 


0.022 


1.666 


0.022 


Abo 


17 
Cj 


4.623 


0.251 


1.922 


4.084 


A IK) 


17 


O/IO 

0.z4z 


U.oUz 


1 1 fin 


U.oUO 


A72 


E 


4.972 


0.013 


2.517 


0.044 


A '7c;a. 
A75 


17 
Cj 


5.463 


0.301 


2.580 


0.304 


A TT 


17 


A Q1 c; 
4.0IO 


U.U04 


0. ilD 




A '70 

A7o 


17 

Cj 


4.837 


0.015 


3.205 


0.051 


A QA 

Ao4 


17 

XL* 


mo 


U.UDl 


1 ani 


U.140 


AnomzUrL 


17 
Cj 


4.360 


0.205 


1.903 


3.410 


rJU+oUoo ' 


17 

Ej 


4.705 


0.113 


3.133 


0.090 




17 
Cj 


4.394 


0.213 


3.921 


3.546 


■p»^ 1 a 
JJal 


17 

Cj 


5.189 


0.301 


1.295 


0.302 


JJerltz 


T7 

Cj 


4.369 


0.011 


1.831 


0.044 


T7/^T3 /id 


17 


4.616 


0.023 


2.278 


0.031 


HaTrl4^ 


E 


4.225 


0.164 


2.473 


2.695 


Hel-2d 


E 


4.475 


0.038 


3.867 


0.097 


He2-429^ 


E 


4.362 


0.205 


3.771 


3.418 


He2-430^ 


E 


4.378 


0.209 


4.054 


3.485 


He2-447^ 


E 


5.031 


0.301 


2.773 


0.301 


Hul-l 


E 


5.055 


0.066 


2.686 


0.176 


IC 1295 


E 


4.993 


0.301 


1.817 


0.301 


IC 2003 


E 


4.954 


0.028 


3.445 


0.109 


IC 2149^'*! 


E 


4.483 


0.008 


3.014 


0.023 


IC 351 


E 


5.127 


0.114 


3.066 


0.218 


J320'^ 


E 


4.772 


0.012 


3.700 


0.074 


Jnl 


E 


4.980 


0.301 


1.844 


0.301 


JnErl 


E 


5.010 


0.010 


1.191 


0.034 


Kl-16 


E 


5.690 


0.264 


1.786 


0.298 


Kl-17^ 


E 


5.343 


0.301 


1.517 


0.302 


K2-l^ 


E 


5.954 


0.301 


0.663 


0.301 


K2-5^ 


E 


4.813 


0.034 


1.611 


0.015 


K3-5^ 


E 


5.375 


0.301 


2.644 


0.302 



Table 2 — Continued 



name morph. log Teff AiogT log L/Lq AiogL 



K3-21'* 


E 


4.392 


K3-57^ 


E 


5.441 


-ivO-Ul 


E 


4 299 


K3-64^ 


E 


4.301 


K3-68'' 


E 


5.476 


K3-76^ 




4.265 


K3-79^ 




4.230 


K3-80* 




4.272 


K3-82'' 


E 


4.427 


K3-92 


E 


4.901 


K4-5* 


E 


5.329 


Ml-2 


E 


4.736 


Ml-6 




4.780 






4.500 


Ml-66^ 


E 


5.037 


Ml- 73 


E 


4.744 


Ml-78* 


E 


4.375 


M2-2* 


E 


5.216 


M2-40* 


E 


4.418 


M2-44^ 


E 


5.307 


M2-45°'''' 


E 


4.279 


M2-50^ 


E 


4.975 


M3-3^ 


E 


5.022 




E 


4.367 


M4-ll^ 


E 


5.359 


Me 1-1 


E 


4.806 


NOP 4n'^''i 


E 


4.476 




E 






E 


4.711 




E 


5.075 




Hi 


4.0 ( U 


NGC 6210'= 


E 


4.814 


NGC 6543'=''* 


E 


4.690 


NGC 6720 


E 


5.148 


NGC 6741'' 


E 


5.356 


NGC 6781'= 


E 


5.022 


NGC 6826'='<l 


E 


4.511 


NGC 6853 


E 


4.777 


NGC 6905^= 


E 


5.117 


NGC 7008 


E 


5.701 


NGC 7048^ 


E 


5.256 


NGC 7354^= 


E 


4.989 


NGC 7662'= 


E 


5.051 


PB9* 


E 


4.988 


PC19^ 


E 


4.949 


Pel-15'i 


E 


4.422 


Pel-16 


E 


4.886 


Pel-20a 


E 


5.608 



0.049 


2.760 


0.159 


0.301 


3.264 


0.302 


0.188 


3.537 


3.126 


0.189 


2.871 


3.137 


0.301 


2.436 


0.302 


0.178 


4.240 


2.939 


0.166 


3.590 


2.728 


U. loU 


o.oiU 


z.yoU 






o.DO ( 


0.006 


2.737 


0.018 


U.oUi 


i. i\JO 


U.oUz 


0.024 


2.384 


0.155 


0.014 


3.273 


0.092 








0.301 


3.080 


0.301 


U.UlU 


o.yoo 




n one 


z.cSoi 




0.301 


3.229 


0.302 


0.218 


3.971 


3.630 


0.301 


3.191 


0.302 


0.182 


4.081 


3.018 


U.oUl 




U.oUi 


0.302 


2.332 


0.306 




z.cSyi 


o.44i 


0.301 


2.153 


0.302 


0.030 


3.969 


0.164 


U.Ui ( 


Q 01 O 


U.U44 


0.020 


3.386 


0.083 


0.018 


4.i7o 


0.062 


U.Uiy 


Q A 

O.40Z 


U.U4y 


0.007 


3.828 


0.034 


0.014 


3.477 


0.086 


0.051 


3.080 


0.045 


0.026 


2.856 


0.059 


0.224 


2.870 


0.153 


0.040 


2.104 


0.112 


0.034 


3.321 


0.074 


0.016 


1.348 


0.102 


0.042 


3.668 


0.093 


0.254 


3.063 


0.301 


0.301 


2.222 


0.304 


0.301 


3.613 


0.301 


0.069 


3.395 


0.184 


0.301 


3.273 


0.301 


0.301 


3.818 


0.301 


0.036 


3.616 


0.111 


0.029 


4.133 


0.136 


0.301 


2.177 


0.303 



Table 2 — Continued 



name morph. log Tefi ^logx log L/L© AiogL 



Pel-21^ 


hi 


5.954 


0.301 


2.610 


0.301 


"Dn 1 


T? 

Hv 


D.vUZ 


U.U4i 


1 Q A1 

l.oUi 


A 1 QA 

U. ioU 


D..od 
rUZ 




4.0U0 




1.0 ( i 


A AQfi 

U.UoD 


Vyi-i 


ill 


4. i oU 


U.UUO 


Q ftQT 
0.00 / 


A C\A 1 
U.U4i 


Vyi-2 


III 


O.U / O 


n r\QA 
U.Uc54 


Q OTK 
O.Z / 


A OQK 
U.ZOO 


weoDo 


III 


A f^OO 

4.oyy 


U.z4o 


1 Q/l /I 
l.O'i-'i 


A AC^A 
4.UOU 


Ait) ' 


TD 
D 


4.ooy 


n 1 


/I A/1 

Z.4U4 


O.OZ ( 


rlU W O 


TD 

D 


4. DO / 


n Aon 


1 oc;o 
l.ZOO 


A AQ1 

U.Uoi 




■R 


4.Z4o 




QQT 


QQQ 

z.ooy 


TT _o yiQ'7a 
rlez-4o r 


rS 


A QTSJ 


n ono 

u.zuy 


o.oio 


o.4cS0 


T^Q 1 a 
ivo-i 


■D 

rS 


a.4oi 


A QA1 




Z.ODO 


A QAO 

U.oUz 


rvo-4 


D 


4.326 


0.196 


2.896 


3.262 


i\o-Z4 


■T) 

rS 


/I OAT 


n 1 c;t 


A 001 


z.ooy 


xVo-v54 


R 

r> 


O.iDZ 




i.oOo 


A QAI 
U.vSUi 


lVo-4o 


rS 


A 01 


n 1 ftA 


O.OOD 


ftO/1 
Z.Dz4 


1VO-4D 


■R 


4.C)44 


n oni 


0011 


'i 'iAA 
0.044 


rvo-ocS 


rS 


/I OQQ 

4.zc5cS 


A 1 

U.lcSo 


Q 1 TO 

o.i * y 


o.Uuo 


ivo-yi 


O 


C OQA 


A QA1 


1 TOA 

i. /yu 


A QAO 

U.oUz 


ivo-yo 


TD 
D 


1 1 Q 


A QA1 


OylO 


A QA1 

U.oUi 


i\.o-y4 


rS 


5.283 


0.301 


2.388 


0.302 


x\.4-0O 


rS 


A Qni 


A 1 QA 

u. ioy 


Q QQC 
0.000 


Q 1 QV 
0. io ( 


iVl 1- o 


■R 


O.ZUO 


U.oUi 


1 1 

z. i 10 


A QAO 
U.oUZ 


iVii-u ( 


R 


c; 1 1 ft 

O.ilD 


A 1 /I T 
U.14/ 


Q TftO 
0. / DZ 


A OT1 
U.Z i 1 


iVii-oy 


R 


o.uyi 


A QA1 

U.oUi 


A AAA 
U.UUU 




A/TI Tea 
Ml- /O 


rS 


5.462 


0.301 


2.154 


0.302 


A/TO /i^a 


R 


A O/lft 


A 1 T1 
U.i / i 


Q /I Ac; 

o.4yo 


QO/1 

z.oZ4 


A/TO /IQa 
iViZ-4o 


R 


A QA1 


A 1 QA 

u. ioy 


o.Oyi 


Q 1 QV 
0. io i 


TV/TO c:oa 


R 


O.DUo 


A QAI 

U.oUi 


ATA 

z.y /U 


A QAQ 

U.oUo 


A/TQ oca 


R 


4.Z0U 


A 1 Tft 

U.i * 


o.ocSi 


ni 

z.yiz 


A/TQ cca 


R 

rS 


/I AfiA 

4.404 


A OOT 


Q Q1 

o.oiz 


Q TftO 

0. (by 


M4-14* 


B 


4.378 


0.209 


3.243 


3.485 


M4-17'' 


B 


5.103 


0.301 


2.125 


0.301 


NGC 650 


B 


5.123 


0.043 


2.162 


0.091 


NGC 2346= 


B 


4.760 


0.011 


3.835 


0.042 


NGC 6778'' 


B 


5.030 


0.070 


2.605 


0.199 


NGC 6881 


B 


4.993 


0.301 


0.000 




NGC 7026"^ 


B 


4.965 


0.029 


2.913 


0.108 


NGC 7027"= 


B 


5.275 


0.301 


5.989 


0.301 


Pcl-14'' 


B 


4.196 


0.152 


4.299 


2.491 


Pcl-17'' 


B 


5.410 


0.301 


2.462 


0.302 


Shl-89=' 


B 


5.221 


0.301 


1.793 


0.302 


Sh2-71^ 


B 


4.890 


0.011 


3.353 


0.041 


A80 


BC 


5.114 


0.014 


1.886 


0.035 


A82'^ 


BC 


4.863 


0.031 


3.071 


0.135 


Hel-1'' 


BC 


5.404 


0.301 


3.203 


0.302 


Hel-6'' 


BC 


5.082 


0.301 


2.102 


0.301 


Hul-2 


BC 


5.058 


0.066 


3.100 


0.174 


K3-36 


BC 


4.890 


0.025 


3.229 


0.120 



-18- 



Table 2 — Continued 



name 


morph. 


log Teff 


AlogT 


log L/L0 


AlogL 


K3-63 


BC 


4.847 


0.301 


3.839 


0.302 




EC 


4.230 


0.166 


3.590 


2.728 




BC 


5.000 


0.301 


2.606 


0.301 


Ml- 7 


BC 


5.045 


0.302 


2.353 


0.307 


Ml-79 


BC 


5.025 


0.302 


2.483 


0.306 


M2-51'= 


BC 


4.776 


0.028 


3.635 


0.153 


M2-53 


BC 


5.035 


0.007 


2.521 


0.023 


M2-55 


BC 


5.031 


0.038 


1.883 


0.099 


NGC 6058 


BC 


4.834 


0.005 


3.435 


0.030 


NGC 6772'' 


BC 


5.147 


0.042 


2.160 


0.094 


NGC 6804'' 


BC 


4.954 


0.007 


3.804 


0.030 


NGC 7139 


BC 


5.020 


0.301 


1.921 


0.301 


PBIO" 


BC 


5.343 


0.301 


3.612 


0.302 



''Temperature, luminosity, and errors for this nebula are from cross-over 
analysis. The temperatures should be considered upper limits. 

''Morphology and dimensions measured from the Schwarz, Corradi, & 
Melnick (1992) images. 

■^Morphology and dimensions measured from the Balick (1987) images. 

■^Temperature, luminosity, and errors for this nebula are from H I 
Zanstra analysis. 



-19- 



Table 3. Physical diagnostics. I. Galactic distribution, extinction, and Zanstra temperatures 



morph. 


<|6|> 


< \A > 


C/3 


< > 


< ^eff.Hell/^'eff.HI > 




[deg.] 


[kpc] 




10^ [K] 




R 


12.6 [12.6] 


0.73 [0.7] 


0.63 [0.73] 


95 [25] 


1.5 [0.6] 


E 


7.2 [7.1] 


0.38 [0.48] 


0.95 [0.77] 


87 [27] 


1.4 [0.6] 


BC 


6.9 [6.9] 


0.37 [0.53] 


1.03 [0.56] 


97 [27] 


1.5 [0.7] 


B 


2.9 [2.5] 


0.21 [0.24] 


1.39 [0.63] 


94 [29] 


1.2 [0.5] 


R + E 


8.9 [9.6] 


0.59 [0.5] 


0.84 [0.77] 


90 [26] 


1.43 [0.6] 


B + BC 


4.1 [6] 


0.21 [0.24] 


1.26 [0.63] 


96 [27] 


1.3 [0.6] 



^AlogTeff < 0.1, He II Zanstra temperatures only. 
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Table 4. Physical diagnostics. II. Abundances 



morph. <He/H>« < 0/H > < N/H > < Ne/H > < Ar/H > 



R 1.10 (18) 3.61 (16) 1.71 (14) 0.96 (16) 1.29 (7) 

E 1.20 (33) 3.34 (32) 1.51 (31) 0.87 (32) 1.28 (27) 

BC 1.39 (9) 3.22 (8) 1.92 (6) 1.36 (5) 0.60 (4) 

B 1.49 (16) 4.05 (13) 5.66 (13) 0.94 (11) 1.36 (9) 



^He/H abundances are multiplied by 10; 0/H, N/H, and Ne/H abun- 
dances are multiplied by 10^; Ar/H abundances are multiplied by 10^. 
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Fig. 1. — Apparent spatial distribution of the sample PNs (Galactic longitude versus Galactic 
latitude). Open circles: R PNs; asterisks: E PNs; filled squares: B PNs; open squares: BC PNs. 
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Fig. 2. — Logarithmic distribution of the scale height for the different morphological types. Upper 
left: round PNs; lower left: upper right: bipolar PNs; lower right: bipolar core PNs. 
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Fig. 3. — Central stars of Planetary Nebulae in the HR diagram. Symbols as in Figure 1. The 
constraints in the temperature and luminosity errors, in this figure, are ATgff < 0.1 dex and 
AL < 0.3 dex. We use the He ll Zanstra temperatures only. Evolutionary tracks are for M=0.55, 
0.6, 0.7, 0.8 0.9, 1, 1.2, and 1.4 Mq, from Stanghellini & Renzini (2000). 
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Fig. 4. — Mass distribution of the CSs, from the sample of Figure 3. Upper left: CSs of round PNs; 
lower left: elliptical PNs; upper right: bipolar PNs; lower right: bipolar core PNs. 
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Fig. 5. - Same plot than in Figure 3, except with ATeff < 0.4 dex and AL < 0.8 dex. We use 
the He II Zanstra temperature and, when those are not available, the H I Zantsra temperatures. 
Symbols and tracks as in Figure 3. 
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Fig. 6. — Mass distribution of the CSs, from the sample of Figure 5. Upper left: CSs of round PNs; 
lower left: elliptical PNs; upper right: bipolar PNs; lower right: bipolar core PNs. 
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